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Abstract 

The isospin distribution of particles and fragments in collisions Ru, 

and at beam energy 400 AMeV is studied with 

isospin dependent QMD model. We find that the rapidity distribution of differential 
neutron-proton counting in neutron rich nucleus-nucleus collisions at intermediate 
energies is sensitive to the isospin dependent part of nuclear potential. The study 
of the N/Z ratio of nucleons, light charged particles (LCP) and intermediate mass 
fragments (IMF) shows that the isospin dependent part of nuclear potential drives 
IMF to be more isospin symmetric and emitted nucleons to be more neutron rich. 

From the study of the time evolution of the isospin distribution in emitted nucleons, 

LCP and IMF we find that neutrons diffuse much faster than protons at beginning 
and the final isospin distribution is a result of dynamical balance of symmetry 
potential and Coulomb force under the charge conservation. 

PACS numbers: 25.70.-z, 24.10.Cn 
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Nuclear multifragmentation has been received more and more attention in both theo¬ 
retical and experimental studies since people believe that it carries abundant information 
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of nuclear dynamics and the properties of nuclear matter, especially, the possible associ¬ 
ation with liquid-gas phase transition. Following the establishment of radioactive beam 
facilities at many laboratories, the study of nuclear multifragmentation for isospin asym¬ 
metric nuclear system becomes possible. For isospin asymmetric systems, a description 
of bi-component system should be introduced while a bi-component system will mani¬ 
fest a richer thermodynamical behavior. The new features of phase separation process 
for asymmetric nuclear system have been discussed earlier by Barranco and Buchler [1] 
with phenomenological equation of state and later by Muller and Serot [2] with relativis¬ 
tic mean held model. It was pointed out in [1, 2] that the spinodal is dehned by the 
chemical instability rather than by mechanical instability for isospin asymmetric systems 
, which allows gas and liquid to have different concentration of protons and neutrons. 
The dynamical properties of asymmetric nuclear matter were investigated in [3, 4, 5] and 
it was predicted that the spinodal decomposition is accompanied by a collective diffu¬ 
sion of protons from low- to high-density region. Experimentally, it was found that the 
N/Z ratio of gas signihcantly exceeds that of liquid in the multifragmentation process 
in heavy ion collisions [6, 7, 8]. It was further showed that the properties of liquid-gas 
phase transition depend on the value of the symmetry energy coefficient[2, 5]. As is well 
known, the symmetry energy term of EOS ( both the density dependence of symmetry 
potential energy and the value of symmetry energy coefficient) has large uncertainties. 
For example, the theoretically predicted value of symmetry energy coefficient is about 
27 — 38MeV by non-relativistic Hartree-Fock approach [9], 35 — AQMeV by relativistic 
mean held approach [10, 11, 12], SlMeV by Brueckner-Hartree-Fock theory (BHF) [13] 
and 28.7MeV by extended BHF theory [14]. Furthermore, recent study has showed that 
the symmetry energy coefficient increases as the isospin asymmetry increases and the in¬ 
creasing slope is rather different for different version of Skyrme force [15]. Therefore it 
might be interesting to study the sensitivity of various observables of nuclear multifrag¬ 
mentation to symmetry energy coefficient. We will hrst investigate the dependence of the 
various observables in multifragmentation on the symmetry energy coefficient and try to 
hnd the sensitive observables to the symmetry energy coefficient. Then the effects of the 
isospin dependent part of nuclear potential on the isospin distribution of emitted nucle¬ 
ons, LCP and IMF in multifragmentation are studied. In order to study the dependence 
of the isospin distribution on the isospin asymmetry of the system we perform a set of 
calculations for the mixing reactions of four mass 96-1-96 systems Ru-|-Ru, Zr+Zr, Ru-|-Zr, 
and Zr-|-Ru with corresponding N/Z ratios of 1.18, 1.4, 1.28, 1.28. It is noticed that 
is of 6 neutron dehciency and is of 6 neutron excess compared with the most stable 
isotopes and . 

The isospin dependent quantum molecular dynamics (IQMD) model [16, 17] is used 
in the calculations. In this model, each nucleon is represented by a Gaussian wave packet 
centered at ri and pi. The time evolution of ri and pi is governed by Hamiltonian equation 
of motion: 

dH dH 



The Hamiltonian H consists of the kinetic energy and effective interaction potential 
energy, 

H = T + U. (2) 

The effective interaction potential energy includes the nuclear local interaction potential 
energy and Coulomb interaction potential energy, 

U = Uloc + Ucoul- (3) 


The local interaction potential energy can be obtained by 


Uloc - 



(4) 


where Vioc is the potential energy density which can be obtained from Skyrme type inter¬ 
action and reads as: 


^ a pjrf P pijf Cs {ppjr) - pn{r)f 

2 po 3 pI 2 po 


(5) 


The third term in the right side of equation (5) is the symmetry potential energy 
density, in which a linear density dependence of for the symmetry potential energy is 
adopted . In general, Vgym can be expressed as: 

Vsym — "o pU{u). (6) 

z p 

Here F{u) gives the density dependence of symmetry potential energy and u = p/po . 
Phenomenologically, Prakash, Ainsworth, and Lattimer proposed the forms of F(u) as u, 
and 2u ^/(1 -|- u) [18],respectively. For the form of or 2v ?/(1 -|- u), we can hardly 
get the analytical expression of Usym from expression (4). However, it has been shown 
in many microscopic studies that the linear density dependence of symmetry energy was 
almost valid at not very far from normal densities, (see, for example, [13, 14]). Thus in 
this work a linear density dependence of symmetry potential energy, F{u) = u, is still 
adopted for the convenience of performing IQMD calculations. Of course the knowledge of 
the density dependence of the symmetry potential energy is very much concerned, which 
will be studied in the future. The symmetry potential energy term in equation (5) can be 
re-written as 

V.Up.S)^U(^-)S"p. (7) 

^ Po 

where 


^ _ Pn ~ Pp 
Pn + Pp 

The corresponding symmetry energy coefficient is 

3 / 2 /s \ n Us Cs 

a-sym = g(2 — l)e^ 4- — ~ 4- —. 


( 8 ) 

( 9 ) 
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Here e% is the Fermi energy of the symmetric system at p = po and taken to be 38 MeV. 
The effect of the isospin dependent part of nuclear potential on the isospin distribution in 
multifragmentation is tested by taking different values of in the calculations. In this 
work we take Cs equals 0 MeV, 27MeV, 35MeV and SOMeV, corresponding to the values 
of symmetry energy coefficient agym'- 0 MeV,27MeV, 31MeV and 38MeV,respectively. 

The experimental isospin dependent binary elastic scattering cross section is used [19] 
in calculations. In [19] it was shown that up to hundreds MeV the elastic proton-neutron 
cross section is about 2 — 3 times larger than that of proton-proton (neutron-neutron)’s. 


Concerning the treatment of the Pauli blocking, we hrstly distinguish protons and 
neutrons, and then we use the following two criteria: 


and 


ATT n ATT n 
—> — 
q ' q — /I ’ 


( 10 ) 


Puod, = l-{l- fj), (11) 

where /j is the distribution function in phase space for particle i and reads as 

fi{r, p,t) = ^ exp(-(r - n {t)f/2L^) exp(-(p - ^ {t)f2L^/h^), (12) 

where L is a parameter which represents the spatial spread of wave packet, and r* (t) and 
Pi {t) denote the center of the wave packet of particle i in coordinate and momentum space 
respectively. The hrst condition(expression (10)) gives the criterion for the uncertainty 
relation of the centroids of Gaussion wave packets of two particles. The second one 
is the probability of the Pauli blocking effect for the scattering of two particles, which 
is especially useful for collisions of heavy nuclei. In this paper, we construct clusters in 
terms of the conventional coalescence model [20], in which particles with relative momenta 
smaller than Pq and relative distances smaller than Rq are considered to belong to one 
cluster. In this work Rq and Pq are taken to be 3.5 fm and 300 MeV/c,respectively, 
following Ref. [21]. In addition, only the cluster with reasonable proton number Z and 
neutron number N are selected in order to get rid of nonphysical clusters. 

The soft EOS {K = 200MeV) is used in the calculations, and the corresponding main 
parameters are listed in Table 1. 

First of all, in order to check our IQMD model, we calculate the rapidity distribution of 
protons to make comparison with experimental data of [23]. Fig. 1 shows the calculated 
rapidity distribution of protons (including deuterons for comparing with experiments) 
for four colliding systems: Zr-|-Zr (solid line), Ru-|-Ru (dashed line), Zr-|-Ru (dotted 
line) and Ru-|-Zr (dash-dotted line) aX E = dOOHMeV and b = 0/m ,in which Cs = 
35MeV is used and the angle selection (10° < diab < 28° and 34° < 9iab < 145°) is taken 
which is the same with that of the experimental data [23]. The experimental data are 
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also given in the figure (solid circles and open circles refiected from experimental data 
assuming a backward/forward symmetry). This figure shows that our calculation results 
can reproduce the experimental data quite well, only at the projectile and target region, 
the calculated results deviate from experimental data, which is because the calculated 
results are for head on collisions while the experimental data are for central collisions. It 
is also shown in the figure that the rapidity distribution of protons for Ru+Ru is higher 
than that for Zr+Zr at whole rapidity region because of the difference of their initial 
proton number. The rapidity distribution for asymmetric reactions Zr+Ru and Ru+Zr is 
between that of Zr+Zr and Ru+Ru. It is checked that the results do not change much 
when Cs is taken to be from 27MeV to 50MeV. 


Fig. 1 


Now let us investigate the sensitivity of the proton number counting to the value of Cs 
because the proton number counting for the mixing reactions of four mass 96+96 systems 
Ru+Ru, Zr+Zr, Ru+Zr, Zr+Ru was experimentally measured by FOPI/GSI recently 
[22, 23]. In [23] a normalized proton number counting Rz is introduced: 


Rz 


2*Z - - Z^^ 

^Zr _ ^Ru 


(13) 


By definition Rz = 1 for Zr+Zr, Rz = —1 for Ru+Ru. Fig. 2 shows the rapidity 
distribution of Rz calculated with Cs = OMeV, 27MeV, 35MeV and bOMeV. And 
we can see that Rz can explore the non-equilibrium effect in nucleus-nucleus colliding 
process successfully. However the results with Cs= QMeV, 27MeV, 35MeV and bOMeV, 
respectively, are indistinguishable. It means that the normalized proton number counting 
Rz is not sensitive to the value of symmetry energy coefficient, which may be understood 
from the definition of Rz where the effect of changing on Rz is significantly reduced. 
So we have to seek the other observables which are sensitive to Cs- 


Fig. 2 


Here we introduce the normalized differential neutron-proton counting which is defined 
as: 


< n > — < p > 

< n > + < p > 


(14) 


Where < n > and < p > is the average number of emitted neutrons and protons, respec¬ 
tively. The relation between Nnp and average N/Z ratio of emitted nucleons is 


N/Z 


1 + A'np 
f-Nnp' 


(15) 


The rapidity distribution of normalized differential neutron-proton counting has been 
shown in Figs. 3 with different Cs of f)MeV, 27MeV, 35MeV, and bOMeV, in which a). 
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b), c), d) corresponds to reactions Zr+Zr, Ru+Ru, Zr+Ru, and Ru+Zr, respectively. The 
dependence of N^p on Cs is very pronounced at the projectile and target rapidity region. 
The general tendency is that N^p increases with Cs- When the isospin dependent part 
of nuclear potential is switched off {Cs = 0), the rapidity distribution of becomes 
flat for symmetric reactions Zr+Zr and Ru+Ru while for asymmetric reactions Zr+Ru 
and Ru+Zr, it becomes an inclined line from Zr side to Ru side, which indicates the 
non-equilibrium effect in reactions. And the N/Z ratio of emitted nucleons (by using 
Eq. 13) at central rapidity {Y = 0) is almost equal to the initial N/Z ratio. As soon 
as the isospin dependent part of nuclear potential is switched on Nnp increases with Cs 
for off central rapidity cases. And one can End from Fig. 3 a) and Fig. 3 b) that at 
projectile and target rapidity region the difference in N^p calculated with different Cs is 
more pronounced for Zr+Zr than that for Ru+Ru. It indicates that the sensitivity to 
Cs increases as the system becomes more neutron rich, which can be understood from 
the expression of the symmetry potential energy term(see Eq.(7)). Therefore the rapidity 
distribution of normalized differential neutron-proton counting in neutron rich nuclear 
collisions at intermediate energies can be considered to be a good candidate of sensitive 
observables to the symmetry energy coefficient. 

Fig. 3 

The N/Z ratio of particles and fragments is relevant to the mechanism of multifrag¬ 
mentation for isospin asymmetric nuclear systems. So we further study the N/Z ratio of 
emitted nucleons, LCP (the cluster of Z=1 or 2) and IMF (the cluster of Z=3-16) as well as 
the average N/Z ratio in regardless of fragments or particles at projectile (|R — 1| < 0.5), 
central (|R| < 0.5) and target (|R + 1| < 0.5) rapidity regions, respectively. Fig. 4 shows 
the calculated results with different Cs for reaction Zr+Zr. The general feature of the 
N/Z ratio of particles and fragments is: the largest one is the N/Z ratio of nucleons, then 
that of LCP, and smallest is that of IMF among them. One can further find from the 
figure that the symmetry potential drives the N/Z ratio of LCP and IMF to approach to 
unit from sides of N/Z > 1 and N/Z < 1, and simultaneously the N/Z ratio of nucleons 
is driven to the higher value than initial one, i.e., to be more neutron rich. The larger the 
Cs is the stronger the effect is. It implies that the diffusion of neutrons increases faster 
than that of protons as Cs increases. However, the total number of neutrons and protons 
has to be conserved and there is a balance of the numbers of neutrons and protons in 
nucleons, LCP and IMF at central, projectile and target rapidity region. Therefore the 
N/Z ratios of all products in multifragmentation is interdependent. If one looks at Fig. 4 
more carefully one can find that at the central rapidity region, the N/Z ratio of LCP at 
central rapidity region decreases while the N/Z ratio of nucleons increases as Cs increases 
but the former one is more sensitive to Cs than the later one. On the other hand, the N/Z 
ratio of IMF increases towards unit when the isospin dependent term of nuclear potential 
is switched on. But generally, the dependence of N/Z ratio of IMF on Cs at projectile 
and target rapidity region is a little weaker than that of nucleons at projectile and target 
region and LCP at central rapidity region. Therefore we may also consider the N/Z ratio 
of LCP at central rapidity region to be a relatively sensitive observable to Cs- The similar 
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tendency of the sensitivity of N/Z ratio of emitted nucleons, LCP and IMF to Cs is found 
for reaction Ru+Ru. 


Fig. 4 


It is convenient for a transport model to study the time evolution of isospin distribu¬ 
tion, from which we can also understand more deeply the dynamical balance of isospin 
distribution in emitted nucleons, LCP and fragments. In £g. 5, we show the time evolu¬ 
tion of the average N/Z ratio of all particles and fragments a), nucleons b), LCP c), and 
IMF d). The time ranges from 30/m/c to 100/m/c. Here 30/m/c is the time when the 
density in a sphere of Ifm^ at center drops to po- In Fig. 5, the corresponding N/Z ratios 
at different time are given by different line types. From this figure, we hrstly see that 
at 30/m/c the N/Z ratio for nucleons and LCP is larger than unit, especially the former 
one. But for IMF it is much smaller than unit( it should be noticed that the number of 
IMF at time is small). It means that neutrons diffuse much faster than protons, which 
simultaneously results in emitted IMF being neutron deficient . As time goes on, the 
N/Z ratio of nucleons goes down which means that more protons are emitted during this 
period. After 50/m/c, the N/Z ratio of emitted nucleons decreases very slowly. The final 
N/Z ratio of emitted nucleons at central rapidity zone closes to the initial N/Z ratio of 
colliding systems but still higher than initial one. While at projectile and target rapidity 
region it is still much higher than initial N/Z ratio. And at the same time span, the 
N/Z ratio of LCP and IMF firstly goes up far away from unit before 50/m/c because of 
emitting of more protons during 30/m/c - 50/m/c, and then goes downward to approach 
to unit. So from the above discussion about the time evolution of the isospin distribution 
in multifragmentation for heavy ion collisions at the energy studied we can draw two con¬ 
clusions: 1) neutrons diffuse much faster than protons so the number of neutrons emitted 
much larger than that of protons at the beginning , and 2) the final isospin distribution is 
a result of dynamical balance of symmetry potential and Coulomb force under the total 
charge number conservation which leads to the N/Z ratio of IMF more close unit and 
emitted nucleons more neutron rich. 


Fig. 5 


In summary, in this work we have introduced a normalized differential neutron-proton 
counting for rich nuclei heavy ion in collisions at intermediate energies as a sensitive 
observable to probe the symmetry energy coefficient. We have also studied the dependence 
of the average N/Z ratio of nucleons, LCP and IMF at central, projectile and target 
rapidity region on (75 and we find that at the energy studied in this work the general 
effect of isospin dependent part of nuclear potential is to drive emitted nucleons to be 
more neutron rich and IMF to be more isospin symmetry. The stronger the isospin 
dependent part of nuclear potential is the stronger the effect is. It is also shown that 
the N/Z ratio of LCP at central rapidity region is relatively more sensitive to Cs than 
the others. The time evolution of the isospin distribution in multifragmentation process 
shows that neutrons diffuse faster than protons and the final isospin distribution is a result 
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of dynamical balance of symmetry potential and Coulomb force under the total charge 
number conservation. 

Table. 1 Parameters used in calculations 


a{MeV) 

l3{MeV) 

7 

Poifm ■^) 

K{MeV) 

L{fm) 

CyukiMeV) 

-356 

303 

7./6. 

0.168 

200 

2.1 

-5.5 
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CAPTIONS 


Fig. 1 The rapidity distribution of protons and deuterons for Ru, Zr, 

Ru, ^^Zr Zr at E = dOOAMeP, b = Ofm. The experimental data of 
[23] are also given in the hgure. 

Fig. 2 The rapidity distribution of proton counting number Rz with different Cs for the 
same reactions given in Fig. 1. 

Fig. 3 The rapidity distribution of differential neutron-proton counting with different 
Cs for a) -F®® Zr, b) ^^Ru -F®® Ru, c) ^^Zr -F®® Ru, and d) ^^Ru -F®® Zr at 400 
AMeV, b = Ofm. 

Fig. 4 The average N/Z ratio of emitted nucleons, TCP and IMF at projectile rapid¬ 
ity region, central rapidity region, and target rapidity region with different Cs for 
^^Zr -F®® Zr at 400 AMeV, b = Ofm. 

Fig. 5 The time evolution of average N/Z ratio of emitted nucleons, TCP and IMF at 
projectile rapidity region, center rapidity region, and target rapidity region for the 
same reactions as Fig. 4. 
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